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We study the phase diagram of the Hubbard model in the weak-coupling limit for coexisting 
spin-density-wave order and spin-fluctuation-mediated superconductivity. Both longitudinal and 
transverse spin fluctuations contribute significantly to the effective interaction potential, which cre¬ 
ates Cooper pairs of the quasi-particles of the antiferromagnetic metallic state. We find a dominant 
d x 2 _ y 2 -wave solution in both electron- and hole-doped cases. In the quasi-spin triplet channel, the 
longitudinal fluctuations give rise to an effective attraction supporting a p-wave gap, but are over¬ 
come by repulsive contributions from the transverse fluctuations which disfavor p-wave pairing com¬ 
pared to d x 2 _ y i. The sub-leading pair instability is found to be in the g- wave channel, but complex 
admixtures of d and g are not energetically favored since their nodal structures coincide. Inclusion of 
interband pairing, in which each fermion in the Cooper pair belongs to a different spin-density-wave 
band, is considered for a range of electron dopings in the regime of well-developed magnetic order. 
We demonstrate that these interband pairing gaps, which are non-zero in the magnetic state, must 
have the same parity under inversion as the normal intraband gaps. The self-consistent solution to 
the full system of five coupled gap equations give intraband and interband pairing gaps of d x 2_ y 2 
structure and similar gap magnitude. In conclusion, the d x 2_ y 2 gap dominates for both hole and 
electron doping inside the spin-density-wave phase. 

PACS numbers: 74.72.-h,75.25.-j,75.40.Gb,78.70.Nx 


I. INTRODUCTION 

Pairing of electrons by exchange of spin fluctua¬ 
tions is a popular paradigm proposed for unconventional 
superconductivity including Fe-based superconductors, 
cuprates, and heavy fermion systems. Since many of 
these systems exhibit an ordered magnetic phase coex¬ 
isting with, or in close proximity to, the superconducting 
phase, a small number of studies have addressed the sub¬ 
sidiary problem of pairing of quasi-particles in the sym¬ 
metry broken spin-density-wave (SDW) phase. While 
this problem has a long history, recent experimental and 
theoretical developments have led a number of authors 
to revisit it' 1 3 

For the simplest case of a doped one-band Hubbard 
model with standard commensurate (ir, it) ordering, the 
existence of pairing and its consequences for the symme¬ 
try of the superconducting order parameter was initially 
investigated in a seminal paper by Schrieffer et al.f 4 where 
the effective pairing interaction was obtained within the 
random phase approximation (RPA) arising from lon¬ 
gitudinal spin fluctuations in the magnetically ordered 
phase. These authors neglected the contribution to pair¬ 
ing from the transverse spin fluctuations corresponding to 
the Goldstone mode of the spin symmetry broken state, 
arguing that while such modes lead to a divergent contri¬ 
bution to the spin susceptibility at the ordering vector Q, 
the coherence factors of the SDW phase screen the bare 
electron-electron interaction vertex, which therefore van¬ 
ishes at Q. Soon after, Frenkel and Hanke 5 showed that 
the transverse fluctuations do contribute to the pairing 
interaction in the same order as the longitudinal fluctua¬ 


tions; the divergence of the transverse spin susceptibility 
is eliminated by the coherence factors, but a residual con¬ 
stant interaction remains. 

For electron-doped cuprates, the one-band Hubbard 
model seems to provide a reasonable minimal model 
since the doped electrons reside primarily on the cop¬ 
per sites. Furthermore, calculations of band parameters 
for the electron-doped cuprates point to the fact that 
the Coulomb interaction is smaller than the bandwidth,^ 
in contrast to their hole-doped counterparts. In conse¬ 
quence, the mean field treatment of the SDW order works 
quite well for the normal state properties of the electron- 
doped systems. In particular, ARPES reports show a 
Fermi surface evolution upon increased electron doping 
which agrees well with the band reconstruction of the 
one-band Hubbard model due to commensurate (n, w) 
order. Upon electron doping, the Fermi surface consists 
of electron pockets around (7r, 0), (0,7r). Close to crit¬ 
ical electron doping, the emergence of hole pockets at 
(or close to) the Fermi level around (f, f), (—§, §) may 
occurP 

Theoretically, the study of spin-fluctuation pairing for 
the electron-doped cuprates has been addressed mostly 
in the paramagnetic phased® 3 where the mechanism 
was found to give rise to a gap with d x 2_ y 2 symmetry 
with strong non-monotonic features as a function of mo¬ 
mentum as a result of Fermi surface intersections with 
the magnetic zone boundary, the so-called hot spots. 
At these positions the pairing becomes particularly pro¬ 
nounced. This behavior was found to agree qualitatively 
with ARPESP-' and RarnarP^ observations. In Ref. 22] 
the possibility for a cross-over to d xy symmetry was 


2 


found at large dopings. In a later work, Krotkov and 
Chubukov 16 * studied the spin-mediated pairing gap close 
to the quantum critical point of antiferromagnetic (AF) 
order and found an anisotropic d x 2_ y 2 gap symmetry, 
but the anisotropic behavior was not related to the hot 
spots, as opposed to previous work. A few studies also 
addressed the coexistence of superconductivity and long- 
range AF order but treated the pairing phenomenologi¬ 
cally, and also found a d x 2_ y 2-wave solution. 1 1 tUHl 

For hole-doped cuprates, stronger interactions imply 
that application of a weak-coupling approach to pairing 
is somewhat less justified. In addition, calculations of 
the spin wave spectrum suggest that the commensurate 
( 7 r, 7 r) order in this case is not the ground state solu¬ 
tion)SEfll complicating the theoretical modelling. Spin- 
fluctuation-mediated pairing has, however, been exten¬ 
sively applied to study also hole-doped cuprates within 
the one-band Hubbard model. First, for the paramag¬ 
netic phase Scalapino et al . 2 ^ generalised the approach of 
Berk and Schrieffe/ 22 ! and found a dominant d x 2 ^2-wave 
pairing instability. Later, the study wa s exte nded by var¬ 
ious methods and numerical approached 22 HHI including 
discussions of the possibility of other superconducting 
pairing symmetries arising from spin fluctuations. From 
a strong coupling approach, as derived from the t — J 
model, analysis of the effective pairing at low hole dop¬ 
ing 28 ! suggests a d -wave superconducting ground state in 
any coexistence phase. This is consistent with rigorous 
perturbative weak-coupling calculations in the presence 
of weak density-wave order*® as well as a study of hole- 
doped cuprates in the coexistence phased In addition, 
a recent study of spin-fluctuation mediated superconduc¬ 
tivity in the SDW ordered metal attacked the problem by 
analytical RPA calculations in the large magnetization 
(small pocket) limit 3 and found that interactions were 
dominated by longitudinal fluctuations on the electron- 
doped side, supporting a nodeless d x 2_ y 2 gap, whereas 
on the hole-doped side both longitudinal and transverse 
fluctuations support a nodal d x 2_ y 2 gap. 

Recently, Lu et a/P-Q studied the case of underdoped 
cuprates in a t — J like model and came to a rather dif¬ 
ferent conclusion. In this paper, the coexistence of com¬ 
mensurate AF and superconductivity was investigated 
in a phenomenological model where the pairing interac¬ 
tion arises based on nearest-neighbor magnetic exchange 
neglecting the double occupancy constraint. In the co¬ 
existence phase, this study found the leading supercon¬ 
ducting instability to be triplet p-wave in the case of 
hole doping, providing a potential explanation for recent 
photoemission measurements indicating a “nodal gap”, 
a state with a full gap near the usual positions of the 
d x 2_ y 2 gap no des w hose existence is well-established at 
higher dopings!- 32 * 33 ! It is worth noting, however, that in 
the original strong-coupling study of unconventional su¬ 
perconductivity driven by the spin waves, studied within 
the t — J model, the d x 2_ y 2 -wave symmetry of the su¬ 
perconducting gap was found to be the only stable so- 
lution.^IlSl We also note that there exist other potential 


explanations of the existence of the nodeless gap in the 
literature 

This controversy suggests the need for a better un¬ 
derstanding of the phase diagram of the single-band 
Hubbard model within a single, reliable approximation 
scheme which can encompass paramagnetic, supercon¬ 
ducting and coexistence phases, and which is capable of 
identifying the strength of pairing by both longitudinal 
and transverse spin fluctuations and charge fluctuations, 
and their relative importance for pairing across a large 
doping range, for different electronic structures, and for 
both weak and strong magnetism. Therefore, we extend 
the work of Ref. [3j by a more complete, fully numerical 
solution of the problem to both confirm the analytical 
calculations and extend them to the larger phase dia¬ 
gram. 

We address the question of how the pairing interactions 
and resulting gap symmetry in the one-band Hubbard 
model, treated within the full spin-fluctuation approach 
with self-consistently determined SDW order, evolve as 
a function of doping throughout the phase diagram, in¬ 
cluding the coexistence dome of SDW order and super¬ 
conductivity. We show that the d x 2_ y 2 solution is in fact 
the leading superconducting instability for all electron 
doping levels. In the hole-doped case, we limit ourselves 
to small hole doping only and force the (7r, 7 r) order to be 
stable by a suppression of additional intrapocket nesting 
contributions to the transverse spin susceptibility. This 
approach also yields a coexistence phase with d x 2_ y 2 or¬ 
der. Thus, the weak-coupling approach where pairing is 
mediated by spin fluctuations gives qualitatively different 
results than the ’’strong-coupling” approach of Ref.I 3 T 1 in 
the case of hole doping. We find that the sub-leading in¬ 
stability is a singlet g -wave solution, which shares com¬ 
mon nodes with the d x 2_ y 2 solution along the zone di¬ 
agonals, and has additional nodes along the momentum 
axes. In the hole-doped case where the Fermi pockets are 
located far away from the zone axis, g and d x 2_ y 2 become 
nearly degenerate. However, since the two solutions share 
common nodes, there is no gain in condensation energy 
by e.g. a time-reversal-symmetry broken solution of the 
form d + ig , and the solution of the full gap equation in 
fact favors d x 2_ y 2 over g. 

We find further that extended s-wave is always sup¬ 
pressed in the coexistence phase, in agreement with 
the findings in Ref. [ 3 ] for small doping levels, and on 
the electron-doped side we do not encounter a lead¬ 
ing triplet gap at any moderate doping away from half- 
filling. The spin-fluctuation-pairing mechanism becomes 
strongly suppressed above the critical doping for which 
long-range magnetic order vanishes, since the nesting 
conditions are rapidly weakening as the Fermi surface 
segments of the paramagnetic phase move apart. 

In a recent study of the iron-based superconductors by 
Hinojosa et a/.® it was pointed out that additional inter¬ 
band gaps may develop in the coexistence phase due to 
pairing of two fermions residing on different bands. Such 
pairs are naively expected to be negligible because they 
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involve fermionic states far from the Fermi energy and do 
not alone manifest a Cooper instability. In this regard, 
these interband pairs can be considered as ” anomalous”. 
However, they are non-zero in the SDW state due to a 
combined effect of the SDW gap and coupling to the in¬ 
traband gaps (dubbed in the following as normal gaps). 
In the iron-based system, this effect is predicted to cre¬ 
ate a different superconducting phase which explicitly 
breaks time-reversal symmetry by development of chi¬ 
ral gaps. We find from a self-consistent treatment of the 
coupled gap equations that the one-band Hubbard model 
also supports substantial anomalous interband pairings. 
The structural factor of the anomalous gaps is d x 2_ y 2 as 
for the intraband gaps, but the phase structure of the 
intraband and interband gaps can be different. However, 
unlike the case of iron-based system.®] the chiral solution 
+ e l ^d™2^„2 does not appear to be favored in our 

x y * x y 

numerical calculations. 


II. MODEL AND METHOD 
A. SDW mean field Hamiltonian 


The model is the one-band Hubbard Hamiltonian 
H ^ T U ^ ^ ^ ^ TLim ( 1 ) 

i : j,cr i i,cr 

where c\ a creates an electron on site i with spin a. The 
interaction term U denotes the energy cost associated 
with having two electrons on the same site. In reciprocal 
space the Hamiltonian reads 

H = ^ fikcj^Ckcr + ^ ^ C k'cr C -k'+qCT :C - k +q 3lCko '’ 

kcr k.k ' ,q cr 

( 2 ) 

with 


£k = — 2 t[cos(k x ) + cos(ky)] — At 1 cos(k x ) cos(k y ) — /:t. ( 3 ) 

The parameter —t is the energy gain corresponding to 
hopping between neighboring sites, and — t' denotes the 
energy gain by hopping to next-nearest neighbor sites. 
The doping level of the system is controlled by changing 
the chemical potential fi. The interaction between two 
electrons is first treated in the Hartree-Fock approxima¬ 
tion giving rise to AF ordering of the spins. We therefore 
consider the mean field Hamiltonian 


#sdw = 5Z^ k <r 

k cr 


( ek 

aW \ 

f Qccr A 

l aW 

e k+Q ) 

V ^k+Qcr / 


( 4 ) 


where W = EkK c lc+Qt Ck t) “ ( c k+q; Ck t)] is the 
AF order parameter, and ^ refers to summation over 
the reduced Brillouin zone only. Diagonalization of the 
mean field Ha miltonian lea ds to the energy spectrum 

K ' 0 = 4 ± \J (e^T) 2 + W 2 , 4 = The mag¬ 

netic gap equation is solved self-consistently given the 
hopping integrals t = 1 , t r , the Coulomb repulsion U and 
the doping. A gapless Goldstone mode at Q = (7r,7r) 
in the transverse spin channel is obtained automatically, 
since the magnetic gap equation is equivalent to the con¬ 
dition 1 — URexo _ (Q) = OP 


B. Superconducting pairing interactions 


Higher order interactions in U generate superconduc¬ 
tivity on top of the AF order through longitudinal and 
transverse spin fluctuations, following the original pro¬ 
posals of Refs. [H and EE Since U scatters the bare elec¬ 
trons, the diagrammatics are performed in terms of the 
bare electron Green’s functions, but the Cooper pairing 
takes place between the quasi-particles of the AF state. 
The quasi-particle operators are related to the bare elec¬ 
tron operators by the transformation: 


Dc cr ^k^fccr T ^k/^k(7 7 (^) 

Ck+Qo- = sign(cr)[u k Q!k<T - Uk/Ikcr]- ( 6 ) 

Transverse and longitudinal spin fluctuations give rise 
to fundamentally different interactions. Inspection of 
the interaction vertex formulated in real spac^l shows 
that the charge and longitudinal interaction vertices give 
rise to no spin flips, whereas the transverse interaction 
does. In the latter channel, we have the gapless Gold- 
stone mode of the AF phase which gives rise to a di¬ 
vergent interaction potential between the bare electrons. 
However, when we consider pairing between the quasi¬ 
particles of the AF state, this divergence is removed by 
the coherence factors as noted in earlier worksP^ In the 
case of pairing between opposite spin electrons, spin flip 
processes are possible and the effective interaction is me¬ 
diated both by longitudinal and transverse spin fluctua¬ 
tions. If pairing occurs between same spin electrons, only 
longitudinal spin fluctuations contribute and the pairing 
potential does not include the bare Coulomb repulsion 
U since this acts only between opposite spin electrons. 
The interaction Hamiltonian is formulated in terms of 
the SDW quasi-particles, and in line with earlier work- 4 
we show the interactions in the longitudinal and trans¬ 
verse channel individually, with the transverse part of the 
interaction stated as a spin flip vertex explicitly, 
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H c /z — A Ar Tk k /(aj c , 0 .a^ k ,-a_k ( TQ;kcr + /^k'o-^-k'CT^-ktT/^ko-) + f k k /(a k/o .a^ k ,-/3_kCT/3k<T + /^k' CT ^-k'3 :Q; -ker a k ( T) 


4AT 


kk'cr 


. ( 7 ) 

H± = —TjJTf Tk.k'( a k'cr Q; -k'CT Q; -ko-akg r + /^kV^-k'S^-ktrAc a) + ? k k , (« k , ff a^ k/ 5 I / 3 _kerA cct + / 3 k'o-/^-k'g ;Q! -ka- a kgO > 


kk 'a 


( 8 ) 


Hc/z ~ 9 ^ r k,k'(«k'<T a -k'<T a -k<7ak<7 + ^k'o-^-k'cr/^-kcrAccr) + f k k ' {oi^ a Ot*_w a ft-]s.cr P\lct + /^k' a A-k'cr Q! -ko' a kCT) > 


kk'cr 


(9) 


with 


r k k / = [2C/ - y c (k - k')]Z 2 (k, k') - [2 u - 14 (k - k' + Q)]m 2 (k, k') + V z (k - k')Z 2 (k, k') - V z (k - k' + Q)m 2 (k, k'), 

( 10 ) 

r tk' = V +-(k - k')n 2 (k, k') - R+_(k - k' + Q)p 2 (k, k'), (11) 


r s k,k' = -Vc{ k - k')/ 2 (k, k') - V c (k - k' + Q)m 2 (k, k') - 


The main ingredients in the pairing interactions are 
the spin and charge susceptibilities within the RPA ap¬ 
proximation 


K(q) = 


u 2 xo(<i) 

i + t7xo(q)’ 


(13) 


y*(q) 


u 2 Xo(<i) 

i - Ux o(q)’ 


(14) 


v+- (q) 


u 2 xo (q) 

i-t/ x +-( q )’ 


(15) 


where the spin susceptibilities are defined by 
Xo(q> w ) = 2ff/ &i “‘( TS qW S -c 


and Xo 


jwf dte 1 “ t (TSq(t')Sl q ). As seen from Eqs. (10)- 
( |12[ ) the bare interaction vertices are modified 
by coherence factors of the SDW phase given by 

with 


w 2 (k,k') = §(l + (-!)* 


y == 

VK) 2 +W 2 ,J(e-,) 2 +W2 


) 


u fi = m2 i l 2 iP 2 i n 2 and f/j = (—1,1,1, —1). 

The interaction Hamiltonians stated in Eqs. 0-0 are 
restricted to Cooper pairing between quasi-particles re¬ 
siding in the same band, i.e. (a k(T a:l k and (/3 kCT /?! k3 ;)- 
In section III El we introduce an extended model which 
includes Cooper pairing between fermions residing on dif¬ 
ferent pockets, i.e. mean fields of the form (a k(T /?l k -). 
For now we restrict ourselves to normal intraband Cooper 
pairs and include pair scattering interactions within each 
band as well as between the bands. The expression for 
the interband couplings between pairs of (k,k') with k 
residing on the a band and k' on the (3 band, which 
are dubbed T kk /, T kk / and T k k / are obtained by inter¬ 
changing the coherence factors p 2 (k, k') •<-> Z 2 (k, k') and 


V z (k - k')l 2 (k, k') - 14(k - k' + Q)m 2 (k,k'). (12) 


m 2 (k, k') o n 2 (k, k') in the Eqs. 
with Refs. [3] and 0] 

The presence of magnetic order breaks spin-rotational 
symmetry, but inversion symmetry is preserved. This al¬ 
lows us to express the superconducting gap in an even 
and odd parity form which corresponds to a quasi-spin 
singlet and a quasi-spin triplet gap, respectively. We la¬ 
bel the gaps by the superscript ”s” for even parity (sin¬ 
glet) and superscript ”t” for the odd parity (triplet) gap, 

i.e. A k (!,/t) = (a_i 4 akf) T (a_ktakt)- As we show 
later when discussing the interband (anomalous) pairing 
terms, it is more useful to classify the superconducting 
gaps in the SDW background by parity rather than by 
spin quantum numbers. In this regard calling Cooper 
pairing spin singlet or spin triplet in the SDW back¬ 
ground actually refers to the even or odd parity of the 
wave function, respectively. 

In the derivation of the superconducting gap equations 
the splitting of these two channels leads to the sym- 
metrization/antisymmetrization of the pairing potential 
in the quasi-spin singlet/triplet gap equation 


10 )-(12) in agreement 


A: 


a,(s/t) 


V 

8 N ^ 

k 


, Ir^/b 

[ L k.k' 


,/ti A k ; (s/t) /D“’ (s/t) 

k tanh' k 


D k ; (s/t) 

/8 

, fls/tj ‘-‘k 
' 1 k.k' 


V 

a /3(s/t) 

) — 


- ) 

V 2 knT ) 


l(s/t)Ak; ' ' tanh ( 


2 k B T 

n^ s/t) ■ 


(16) 


and similarly for A k ’^ t ' ) by interchanging a f3. In 
the case of opposite spin interaction the effective pairing 
interactions for the singlet and triplet channel are given 
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by 


rL s / k , - (r k , k ' + 2r+-,) + (r* kik , + 2r+- k ,), (17) 
= (r k , k / - 2r+”) - (r* kik , - 2r+- k< ). (is) 

for the intraband contributions and equivalent ex¬ 
pressions for the interband, T, contributions. Simi¬ 
larly, for same spin electron interactions, the singlet 
and triplet potential is obtained directly by a sym- 
metrization/antisymmetrization of the potential stated 
in Eq. (12 ). 


The strength of the pairings is calculated by evaluating 
the real part of the RPA susceptibilities at zero energy. 
Note that due to broken spin rotation symmetry there 
is a difference between Xo and Xo » The pairing vertex 
V + _(q) diverges at q = Q due to the Goldstone mode 
in the transverse channel. This divergence is removed, 
however, by the coherence factor p 2 (k, k') of the SDW 
phase, as discussed above. 

In the paramagnetic phase where W = 0, pairing takes 
place between the bare electrons and the gap equation 
then reduces to 


A 


(s/t) 


= -ivIX 


(s/t) 


.(s/t) 


AN 


±V' 


' E, 


(s/t), 


with 4 s)/(t) = y^ 2 " 
are given by 


V, 


=u+-[v z (k- k') - 14 (k - k')] ± F+„(k - k'). 


( 20 ) 


C. Gap symmetries in the SDW phase 


In the SDW phase, the effective pairings contain Umk- 
lapp terms, which are the terms in Eqs. (10)-(12) con¬ 
taining the argument k — k' + Q. Due to these terms, 
there are attractive pair scatterings (k, k') on the recon¬ 
structed Fermi surface. This contrasts with the situation 
in the paramagnetic phase, where the spin-fluctuation 
part of singlet potential is always purely repulsive for all 
pairs of (k,k'), as seen from Eq. (201. However, despite 
the presence of a partially attractive potential, a conven¬ 
tional s-wave superconducting gap is not possible. This 
is due to a symmetry constraint on the pairing poten¬ 
tials, which must obey T kjk / = -T k+ Q. k / = -r kk+ Q, 
a property that is fulfilled for both the longitudinal and 
transverse pairing interactions. This symmetry also car¬ 
ries over to the superconducting gap, which must satisfy 


A k = -A 


k+Q- 


( 21 ) 


For the two dimensional square lattice, the gap solutions 
can be classified according to the five irreducible repre¬ 
sentations of the group that are even under reflection 








FIG. 1. (color online) Fermi surfaces in the case of (a) small 
hole doping, (n) = 0.96 and (b) small electron doping (n) = 
1.05 for U = 3 and t' = —0.35. 


through the horizontal plane, i.e. extended s-wave (s*), 
d x 2 _ y 2 , d xy , g and the triplet solution, p x /Py, which is 
doubly degenerate. First we consider whether these solu¬ 
tions all comply with the additional symmetry constraint, 
Eq. (21) present in the SDW phase. This disqualifies the 
d x 


k.k' 

J — TTpr tanh ( 

J p(s/t) V 

2 k B TJ’ 

A\g 



(19) 

Big 




A 2 g 

I 2 - 

The effective 

pairings 

E u 


xy solution and as a result we consider the leading gap 
symmetries from among the set 

cos (k x ) + cos(ky), (22) 

" — ”'~i(k x ) - cos(ky), (23) 

— cos(k v )]sm(k x )sm(k v ), (24) 

P.r. 


s 

d X 2_y2 = cos 

g = [cos(fc x 

sin (k x ), p y 


(25) 


We emphasize that these basis functions are only the low¬ 
est order functions corresponding to the given irreducible 
representations. In the construction of higher order solu¬ 
tions in the SDW phase, it is important to note that an 
odd number of the above basis functions must be multi¬ 
plied in order for the resulting higher order gap function 
to comply with Eq. (211. Further allowed higher order 


solutions can also be achieved by a multiplication of the 
d xy basis function [sin(fe x ) sin^)] with any Ai g , Bi g or 
E u basis functions. A relevant triplet solution is con¬ 
structed by a multiplication of the triplet p-wave with 
the lowest order Ai s and Bi ff basis functions, 

p' x = [cos (k x ) - cos(k y )][cos(k x ) + cos (k y )] sin (k x ), 

(26) 

which for tetragonal symmetry belongs to the same E u 
symmetry representation as the original p-wave solution. 
Therefore, we dub this state p'-wave. Note that in the 
literature on Sr 2 Ru 04 this function is sometimes called 
f x 2 _j, 2 -waveP^ We find that this is the leading solution 
among the triplet solutions for any doping. However, as 
we shall see below, it does not dominate over the singlet 
solutions. 

The consequences of a d x 2_ y 2 solution are manifested 
very differently depending on the Fermi surface geometry. 
In the presence of only electron pockets, as in Fig. Qb), 
it is nodeless at the Fermi surface, as opposed to the 
other solutions, which all display nodes. If both electron 
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and hole pockets are present at the Fermi surface as in 
Fig. E e), the d x 2_ y 2 solution has nodes only at the hole 
pockets, whereas the extended s-wave as well as the g- 
wave solutions exhibit nodes at both types of pockets. 
Naively, we might expect the preferred solution to have 
the minimum number of nodes at the Fermi surface. In 
this respect, the d x 2_ y 2 solution clearly wins, but ulti¬ 
mately the leading solution relies on a detailed investiga¬ 
tion of the structure of the pairing interaction in the SDW 
phase and the solution to the full gap equation. Finally, 
when only hole pockets are present at the Fermi surface 
as shown in Fig. [lja), the p x solution becomes nodeless, 
whereas none of the singlet gap symmetries will provide 
nodeless superconductivity. As we show below, the p x so¬ 
lution is supported from the structure of the longitudinal 
fluctuations, but becomes overall less favorable compared 
to d x 2_ y 2 due to an effective repulsion arising from the 
transverse spin fluctuations. 

D. Linearized gap equation in the SDW phase 


The linearized gap equation does not allow for a 
determination of complex gap solutions, which are 
time-reversal-symmetry-broken (TRSB) solutions, since 
higher order interactions in Ak are removed in the 
linearization process. Therefore, we have addition¬ 
ally solved the full nonlinear gap equation as given in 
Eq. This allows us to compare the real solutions 

and TRSB solutions energetically. The latter type of so¬ 
lutions arise naturally in situations where the pairing po¬ 
tential allows for two degenerate eigenfunctions of AA.k' ■ 
Therefore, they are likely to appear in the triplet chan¬ 
nel, where all solutions to the linearized gap equation are 
two-fold degenerate. Since the linearized equation in fact 
also exhibits accidentally degenerate solutions in the sin¬ 
glet channel upon hole doping, this opens the possibility 
of TRSB states in this case. The prospect of observation 
of transitions to TRSB states as a function of doping or 
other control parameters, potentially the first observation 
of this kind, has recently been the subject of considerable 

attention in Fe-based superconductors 


In order to determine the leading and sub-leading in¬ 
stabilities, we linearize the full gap equation stated in 
Eq. (16) to obtain the eigenvalue problem 


1 


4(2tt) : 


rMk.k'Ak' = AA k q 


(27) 


The eigenvector of this equation 


Al 


A“ 

A k fel 

A ^2 


(28) 


is a function of momentum k. The momentum is lo¬ 
cated either on the electron pocket around (tt, 0), which 
we denote by a, or on the two hole pockets, denoted 
by Phi/Ph. 2 , around (f, f )/(-f, §), respectively, see 
Fig. [lj We solve the eigenvalue problem Eq. (271 for 
k and k' on these three closed pockets, which is equiv¬ 
alent to solving it in the whole magnetic Brillouin zone. 
All intraband and interband interactions are included in 
Mc.k': 


paa 

-pPh 1 a j a 

1 k,k' L k' 

-pfih 2 a ja 

1 k,k 7 L k' 


\ u k' 

|v£> 


-p a ^h 1 i/3 
1 k,k' fc k' 


K%' hl €/K 

j^0h, 2 Ph 1 i(3 /| B I p/3/i 2 Ph 2 j/3 /| B 
L \c' / I u \c' I 1 k,k' t| - / / 1 u ''‘ 


Km'C/ 

-pPhiPh.2 J0 

1 k.k' l k' 


Here, lk> denotes the length of the Fermi surface line seg¬ 
ment represented by the point k' and Ck' is the Fermi ve¬ 
locity at k'. The intraband (Tkk') and interband (Tkk') 
pairings must be expressed in the singlet and triplet sym¬ 
metrized versions in accordance with Eqs. |[l7j-((T8|. The 
largest eigenvalue, A, gives the leading instability since it 
corresponds to the largest superconducting critical tem¬ 
perature and the symmetry of the gap is given by the 
corresponding eigenvector. 


E. Anomalous pairing gaps 

So far, we have restricted the discussion to Cooper 
pairing between quasi-particles residing on the same 
band, which are Cooper pairs of the form and 

(/^ko-Z^-k a)- N° w we introduce an extension of the model 
to also include anomalous pairings between fermions be¬ 
longing to different pockets, i.e. (aj^• The com¬ 
plete interaction Hamiltonian includes pair scattering 
processes between normal and anomalous gaps, as well 
as scatterings between anomalous gaps. 

Anomalous pairs of the form (a involve 
fermions far from the Fermi level, since the two bands 
are gapped by | E k — E k \ > 2W. Nevertheless, the 
anomalous gaps become sizeable due to the coupling to 
the normal intraband pairs and exist only in the SDW 
phase since the coupling between normal and anomalous 
gaps are proportional to W. A detailed analysis of the 
interaction Hamiltonian reveals that even parity intra¬ 
band gaps, i.e. (crk^a^Lk^) — ( a i 4 a -kt) cou pl e to even 
parity anomalous interband gaps. In the iron-pnictide 
study of Ref. 03] the reported anomalous gap was dubbed 
quasi-spin triplet and coupled to an even parity singlet 
intraband gap. However, in that case the triplet gap is 
actually of even parity, which is allowed due to the band 
index. As a matter of fact, in the SDW background it 
is more useful to classify the superconducting gaps by 
parity rather than by spin quantum numbers, as spin- 
rotational symmetry is explicitly broken. Thus, in our 
case as well as in Ref. 03] one has an even-parity normal 
intraband gap coupling to even-parity anomalous pair¬ 
ing contributions. As we shall see below, calculations of 
the normal intraband gaps reveal an even-parity d x 2 _ y 2 - 
wave solution. Because of this, and also as a result of 
the structure of the interaction Hamiltonian, we restrict 
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ourselves to the even-parity channel which we label by a quasi-spin singlet index, s. In the even-parity channel, 

the mean-field interaction Hamiltonian takes the form 


Ha = [a® a (k)a_ k .|.a k f + A^ /3 (k)/3_ k |^ k | + A* /3 (k)[a_ k |/3 k ; - + /J-iq-aiq, - /3_ k ^a k ^] 


+ A «m( k )a-k|4t - A* m (k)/3- ki a kt + h.c. 


(29) 


with the mean fields 

AL(k) = -^E'+ f s k ,k'<& 7 0!w>] + K^k,k'< a L^- k '7 + 4'^-k'7>- (30) 


A ^ (k) = E' *^[ r k,k'(4 7 ^ k ^) + r^, k ^ <«L 7 «t k ^>] - K* r M('>k'^-k' 7 + ^“U 7 )» ( 31 ) 

k' 

A a/3(k) = ^^k^k' ( a k'7 a -k'.5 - ^k^-k'a)’ 

k' 

_ _ 1 r a/3a/3,l,s/ t ot 


(32) 

(33) 

(34) 


A; m (k> = - w Y r^v^W+4»- 2 -(4,« , - k . i >, 

k' 

AUr(k) = - j^E' r ; s k "'’’ 1 ''(4.*4 k .t) + r2S?' ,A ’(4 4 “b't>- 

k' 

In these gap equations we have introduced three new spin-fluctuation-mediated pairing interactions 

K% = T[Vi„(k-k'Hk,k')+y io (k-k' + Q)/x(k, k')] - [ 2 y+_(k-k>(k,k') + 2 U + _(k- k' + Q)i/(k,k')], 

(35) 

K^’ 1 = -[Vio(k - k> 2 (k, k') + V lo {k - k' + Q)n 2 (k, k') + 2U + _(k + k> 2 (k, k') + 2U + _(k + k' + Q)p 2 (k, k')], 

(36) 

r a/3a/3,2 = fc /) + y^ k _ k ' + Q) m 2(^ k ') + 2 U + _ (k + k')m 2 (k, k') + 2U+_(k + k' + Q)Z 2 (k, k')], 


(37) 


where V/ 0 (q) = 2 U — I4(q) + I4(q), and =F in Eq. (35) refers to the quasi-spin singlet and triplet, respectively. In 

- \Oi/3,s -pot/3 . -pa(3 

k,k' — 1 k,k' ' 1 — k,k' 


Eqs. (|30|)-(|34|) the superscript s refers to the symmetrized potentials, e.g. T? k f = r£ k / + T“? This ensures that 


all mean fields are of even parity. 

We have introduced two new coherence factors, which are proportional to the magnetic order parameter W 


W t e k' + e k 




W / 

M (k,k') = _( 


(38) 

(39) 


Note that the coupling between the normal intraband 
gaps and the anomalous pairings occurs due to the pair¬ 
ing stated in Eq. (35) which is proportional to the SDW 
order parameter. 


III. RESULTS 


In the next sections, we discuss the results of spin- 
fluctuation-mediated pairing in the SDW phase. First, 
we restrict the model to normal intraband Cooper pairs 
only. In section |III A| we study how the effective interac¬ 
tions mediated by transverse and longitudinal spin fluc¬ 
tuations evolve as a function of doping. Thereafter, we 
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FIG. 2. (color online) Maximum value of the pairing contri¬ 
butions from longitudinal and transverse spin fluctuations in 
the singlet and triplet channel, r k k ,, in the SDW phase as a 
function of doping, for the parameters t' = —0.35 and U = 3. 


turn to the implications for the evolution of different gap 
symmetry solutions in section IIIB We focus on elec¬ 
tron doping, since this is the most relevant regime for 
a coexistence phase of commensurate AF order and su¬ 
perconductivity. In the case of hole doping, we restrict 
the investigation to small hole doping levels below 10 %. 
Lastly, we provide a discussion of the extended model, 
where anomalous pairings are included, in section uni ] 


A. Pairing interactions from longitudinal and 
transverse fluctuations 

First we focus on the strength and structure of the 
pairing potentials between opposite spin electrons aris¬ 
ing from longitudinal and transverse spin fluctuations 
for hole and electron doping. In Fig. [2| we plot the 
value of the pair scattering T kk , for which the pairing 
strength is maximal for any pair of momenta (k, k') on 
the Fermi surface. In the low doping limit, we observe 
that the interaction through transverse fluctuations is 
much stronger for hole doping than for electron doping. 
For interactions in the longitudinal channel, the situation 
is opposite, since in this case the pairings are strongest 
for small electron doping compared to small hole doping. 

Generally, for the hole-doped system, transverse fluc¬ 
tuations are quantitatively stronger than the longitudi¬ 
nal fluctuations. For electron dopings, the longitudinal 
fluctuations dominate close to half filling, whereas both 
types of pairings contribute more equally at larger doping 
levels. In fact, the strength of the pairing interactions in 
both channels builds up towards the critical electron dop¬ 
ing where the AF order disappears. This naively suggests 
an increase in superconductivity close to critical electron 
doping. However, a closer inspection of the (k, k') de¬ 
pendence of the pairing reveals that it develops a sign 






FIG. 3. (color online) Effective pairing contributions from 
longitudinal (a,c) and transverse (b,d) spin fluctuations in 
the singlet and triplet channel on the hole pocket centered at 
(f, f). The black cross denotes the position of k' and the 

value of r^l, is shown as a function of k around the hole 
pocket. Negative potential contributions correspond to effec¬ 
tive attraction. The zone diagonal is indicated by a dashed 
line. The filling is (n) = 0.96 and U = 3.0, t' = —0.35. Note 
the different color scale for pairing potentials in the longitu¬ 
dinal and transverse channels. 


change for a part of the intrapocket potential and be¬ 
comes strongly repulsive for nearby points at the Fermi 
surface. This leads to an overall gap suppression despite 
the increased interaction strength. 

In general, knowledge of the (k, k') structure of the 
pairing potential is crucial in order to decide what kind 
of superconducting instabilities are favored in the AF 
phase. To address this question, in Figs. [3]and[4]we map 
out both the singlet and triplet pair scattering potentials 
due to longitudinal and transverse spin fluctuations for 
hole- and electron doping, respectively. In the case of 
small hole doping (Fig. [3]), the singlet pairing interac¬ 
tions from longitudinal and transverse spin fluctuations 
are similar in structure; there is an effective attractive 
(negative) interaction for k and k' on the same side of 
the zone diagonal, whereas the interaction becomes repul¬ 
sive (positive) when k and k' are on opposite sides of the 
zone diagonal, as deduced from Fig. [3ja-b). This struc¬ 
ture supports d -wave as well as g -wave states. In this 
case, the pairing interaction due to transverse spin fluc¬ 
tuations is much stronger than the pairing arising from 
longitudinal spin fluctuations, as was also deduced from 
Fig-! and expected from Refs. [51 [51 and flTf, In the triplet 
channel, the effective interaction arising from the two dif¬ 
ferent types of spin fluctuations has very different char¬ 
acter. The longitudinal fluctuations in fact give rise to a 
locally attractive pairing on the entire pocket, as shown 
in Fig. [3](c) . This type of pairing supports a nodeless gap 
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of p -wave character, i.e. Ak cx sin(fc x ). However, due to 
the strong intrapocket repulsive interaction mediated by 
transverse spin fluctuations shown in Fig.0d), the total 
pairing interaction does not allow for a p -wave solution. 3 
Therefore the triplet solution becomes higher order and 
will display nodes at the Fermi level. 

For electron doping, the pairing interaction in the sin¬ 
glet channel is locally attractive for k and k' residing 
on the same electron pocket, as shown in Fig. 0a-b). 
Attractive pairing occurs due to both longitudinal and 
transverse spin fluctuations. As opposed to the hole- 
doped case, it is the longitudinal fluctuations that me¬ 
diate the strongest effective pairing. In this regime, it 
is clearly the d x 2_ y 2 solution which will be favored due 
to the symmetry Tk.k'+Q = —Tk.k'- Upon increasing the 
electron doping the purely attractive intrapocket interac¬ 
tion is replaced by partly repulsive regions and this will 
diminish the resulting superconducting gap value. In the 
triplet channel the pairings from both types of fluctua¬ 
tions will be quite weak, see Fig. 0c-d). As in the case 
of hole doping, the longitudinal fluctuations support a p- 
wave gap since it is attractive for k and k' located at the 
same side of the Fermi pocket and repulsive for k and k' 
on opposite sides, as shown in Fig. 0c). The transverse 
fluctuations display the reverse structure as evident from 
in Fig. 0d), which is the reason why the p-wave solution 
becomes suppressed also on the electron-doped side. 


B. Solutions to the linearized gap equation 


Keeping in mind the structure of the pairing interac¬ 
tions described above, we now turn to the solutions to 
the linearized gap equation Eq. (271 in order to deter¬ 
mine the subleading solutions and the doping evolution 


(a.) T singlet long. 

20 
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(b) T singlet trans. 
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FIG. 4. (color online) Effective pairing contributions from 
longitudinal (a,c) and transverse (b,d) spin fluctuations in the 
singlet and triplet channel on the electron pocket centered at 
(■7T, 0). The black cross denotes the position of k' and the 
value of , is shown as a function of k around the electron 
pocket. The filling is (n) = 1.05 and U = 3.0, t' = —0.35. 
Note the different color scale for singlet and triplet pairing 
potentials. 


of the three leading solutions. We consider three qualita¬ 
tively different types of Fermi surfaces. First, we study 
a Fermi surface consisting of hole pockets at (±f,±^) 
which occurs on the hole-doped side, see Fig. 0a). Sec¬ 
ond, we turn to electron pockets at (±7r, 0) and (0, ±7r) 
which occur for small and intermediate electron doping 
levels, see Fig. 0c), and third, we study the occurrence of 
both electron and hole pockets very close to critical dop¬ 
ing as in Fig. 0e), i.e. the doping level for which AF or¬ 
der disappears. As discussed above, the presence of hole 



FIG. 5. (color online) Solutions to the linearized gap equation 
in the presence of (a,b) only hole pockets (4 % hole doping 
and U = 3), (c,d) only electron pockets (5 % electron doping 
and U = 3) or (e,f) both types of pockets (5% electron doping 
and U = 2.735). The last situation is obtained very close to 
the AF quantum critical point, i.e. for W —> 0. The leading 
solutions are shown in the first column and the sub-leading 
solutions are shown in the second column. In all three cases 
the leading solution is d x 2 _ y 2 , and the sub-leading solution 
is g-wave. The latter differs from s-wave by being odd under 
k x —> —k x and k y —> —k y . Note that in the case of hole dop¬ 
ing the two solutions d x 2_ y 2 and g are nearly degenerate with 
\ d = 0.384 and \ g = 0.381. In the case of small and interme¬ 
diate electron doping, the d x 2_ y 2 solution becomes strongly 
dominant. Upon increased electron doping approaching the 
quantum critical point, the g- wave solution becomes increas¬ 
ingly important, although always subdominant. 
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pockets destabilizes the commensurate AF by additional 
intrapocket contributions to the bare spin susceptibili¬ 
ties. In this case, we force stability of the commensurate 
AF order by turning off the intrapocket contributions by 
hand, which is justified from the work by Chubukov and 
Frenkel 1 ® showing that vacuum renormalization leads to 
a stability of the commensurate AF order at small hole 
dopings. 

By inspection of the intrapocket potential structure on 
the hole pocket shown in Fig. [3^a,b) , we expect either of 
the two solutions d x 2_ y 2 or g for the weakly hole-doped 
system. From the intrapocket potential structure it is not 
possible to qualitatively distinguish between these two so¬ 
lutions. The only other way that the system might choose 
one solution over the other would be from the structure 
of the pairing potential for k on one hole pocket and k' on 
the neighboring pocket. However, from the numerical po¬ 
tential evaluation it turns out that the interpocket pair¬ 
ing contribution is an order of magnitude smaller than 
the intrapocket pairing and its symmetry does not favor 
one of the solutions over the other. This is the reason we 
obtain two nearly degenerate singlet solutions, namely 
d x 2_ y 2 and the g solutions, in Fig. §a,b). 

Turning to the case of 5 % electron doping shown in 
Fig. ,d), we see that the d x 2_ y 2 solution is clearly 
favored over the sub-leading solution g, which becomes 
strongly suppressed. This is a direct consequence of the 
intrapocket attraction between k and k' residing on the 
same electron pocket. Upon increasing doping, the elec¬ 
tron pockets grow in size, but throughout the SDW re¬ 
gion, the d x 2_ y 2 solution continues to be the leading so¬ 
lution although the sub-leading g -wave gets closer. Also 
close to the crossover to the paramagnetic phase, where 
both electron and hole pockets are present at the Fermi 
surface, the leading solution remains d x 2_ y 2 as shown in 
Fig- [5je,f). 


lution is energetically favored. Upon larger hole doping, 
the near-degeneracy of the d x 2_ y 2 and g -wave solutions 
is split and the d x 2_ y 2 solution becomes clearly dominant 
in this regime. In the case of electron doping, the d -wave 
solution is strongest very close to half Filing even though 
this is not where the longitudinal and transverse pairing 
potentials achieve their maximum strengths, see Fig. [2] 
The reason is that in the limit where electron pockets 
are small, the structure of the intrapocket pairing poten¬ 
tials is purely attractive, as shown in Fig.[d|a,b), and this 
strongly supports a d -wave solution. We note that this 
feature of a well-developed gap in the limit of small elec¬ 
tron dopings is an inherent result of the weak-coupling 
approach to the coexistence phase. At critical electron 
doping for which W —> 0, the Fermi arcs just touch the 
magnetic zone boundary and as a consequence, nesting 
by Q on the paramagnetic side is rapidly weakened upon 
increased electron doping. The ordering of the leading 
solutions remains the same as in the SDW phase, with 
the p' solution, which in the paramagnetic phase takes 
the simpler form [cos(fc x ) — cos(fc y )] sin(fca;), the least fa¬ 
vorable. The paramagnetic Fermi surface, which is a 
hole pocket centered at ( 7 r, n) is roughly circular thereby 
preventing nesting not only at Q, but at any q-vector. 
As a result spin-fluctuation-mediated superconductivity 
rapidly dies off. The evolution of the d x 2_ y 2 solution 
was previously discussed in Ref. [SJ where the possibility 
of a different pairing mechanism close to the crossover 
between SDW and paramagnetism was speculated. In 
Ref. [5] the pairing was treated at a phenomenological 
level and here we note that also the full treatment of 
the spin-mediated pairing gives rise to a rapid decrease 
in the superconductivity upon entering the paramagnetic 
region. 

Finally, we show estimations for the critical tempera¬ 
tures of the commensurate AF order, T)v, as well as the 


C. Coexistence phase diagram 


The complete doping evolution of the three leading su¬ 
perconducting order parameters, d x 2 _ y 2 1 g : and p'- wave, 
is shown in the phase diagram in Fig. [6] On the hole- 
doped side, we limit the study to the underdoped region. 
Very close to half filling we observe a near-degeneracy of 
the d x 2_ y 2 and g -wave solutions. The existence of two de¬ 
generate solutions in principle allows for their mixture as 
a TRSB solution, or for orthorhombic distortions as a di¬ 
rect mixture. At the mean field level, the preference of a 
low temperature TRSB solution might be expected when 
it is possible to remove gap nodes at the Fermi surface, 
since removal of gaps from | Ak| leads to a gain in the con¬ 
densation energy. In the present case of the degenerate 
solutions d x 2_ y 2 and g, however, we note that these share 
the same nodes along the zone diagonal. Thus, there will 
be no apparent gain of energy by constructing a TRSB so¬ 
lution of the form d x 2 _ y 2 ±ig. In fact, solving the full gap 
equation as stated in Eq. (161 reveals that the d x 2_ y 2 so- 



(n) 


FIG. 6. (color online) Three leading superconducting instabil¬ 
ities as a function of filling. The SDW region is shown by the 
green area. The three largest eigenvalues to the linearized gap 
equation, Eq. (271, are shown. The d x 2_ y 2 solution (red line) 
dominates at ail fillings. The next-nearest neighbor hopping 
is t = —0.35 and the bare Coulomb interaction is U = 3. 
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FIG. 7. (color online) The critical temperature, T c = 
1.13 e c e 1X , with the energy cut off set to t c = 0.25 and 
t = 400 meV. The superconducting instability is d x 2_ y 2 at 
all fillings. The SDW region which is determined from the 
mean field gap equation is shown by the green area. The 
next-nearest neighbor hopping is t' = —0.35 and the bare 
Coulomb interaction is U = 3. Note that within the weak- 
coupling approach the SDW ground state is always metallic 
away from ( n ) = 1 and therefore, does not cover the Mott in¬ 
sulating behavior of in the underdoped region of hole-doped 
cuprates. 


superconducting ordering temperature, T c in Fig. [7] The 
T c is estimated from the eigenvalue A obtained for the 
leading d -wave solution as a function of doping. Whereas 
T c drops off upon increasing electron doping it shows 
the opposite evolution as a function of hole doping. At 
(n) = 1 the superconducting instability is absent due to 
a full gapping of the Fermi surface by the magnetic order. 
Similar behavior was found by variational cluster pertur¬ 
bation theory in the work by Senechal and coworkersP^ 

As discussed in the Introduction, the results on the 
hole-doped side contradicts the recent findings of Ref. [5TI 
where a p -wave solution appears as the leading instability. 
In fact, the Cooper-pairing in Ref. [31] is mediated by an 
exchange interaction which involves only nearest neigh¬ 
bor sites, while in our approach the interaction potential 
in the Cooper channel invokes effectively also sites which 
are farther apart. The similarity of the treatment of the 
coexistence phase in Ref. [31] and our approach allows for 
a direct comparison of the interaction Hamiltonians. In 
the simplified case of only hole pockets at the Fermi sur¬ 
face, the interaction Hamiltonian reads generally 


Hhoie= r(k,k , )/3L^l k ^_ k ^/3 kV . (40) 

k^k'tr 


In the t -Alike model without double occupancy con¬ 
straint employed in Ref. [31] the effective interaction en¬ 


tering Eq. (40) takes the form 


r t _j(k,k') = - Jtyk 0 k ) [m 2 (k,k') +l 2 (k,k')] 
-J(k + k> 2 (k,k')+p 2 (k,k')], 

(41) 

with J(q) = J[cos(ga;) + cos(< 7 y ] and J > 0. By contrast, 
in the Hubbard model the effective interaction entering 
Eq. (40) is given by 


r H u b (k, k') = re k , ± 2 r+- 


(42) 


with the longitudinal and transverse effective interactions 
stated in Eqs. (101 and The lower sign of Eq. ( [42] ) 
belongs to the triplet channel and this is the source of an 
effective intrapocket repulsion on the hole-doped side as 
shown in Fig. id). On the contrary, the effective interac¬ 
tion as stated in Eq. (41) gives rise to a purely attractive 


triplet potential for k and k' residing on the same hole 
pocket, and therefore the p-wave solution would indeed 
appear to be the dominating instability. Note, however, 
that the model of Ref. [3T] cannot be considered strictly 
as a strong-coupling limit of the single band Hubbard 
model as it does not include a constraint for no double 
occupancies of the fermions explicitly. This could play an 
important role. For example, as mentioned in the Intro¬ 
duction, the original strong-coupling study of unconven¬ 
tional superconductivity driven by the spin waves, stud¬ 
ied within the t-J model with the constraint of no double 
occupancies, does find the da; 2 _. y 2 -wave symmetry of the 
superconducting gap to be the only stable solution! 34 * 35 ! 


D. Anomalous pairing gap 


We introduced the possibility of additional pairing 


gaps in the coexistence phase in Sec. HE Such anoma¬ 
lous pairings were neglected in the previous discussion, 
see Figs. [2][TJ where we demonstrated that the dominat¬ 
ing symmetry of the normal intraband pairing gaps is 
d x 2_ y 2 at all doping values. Now we address the question 
of whether the inclusion of anomalous interband gaps 
introduces essential modifications to this result. The 
anomalous gaps develop under the constraint that they 
must have the same parity as the normal intraband gaps. 
In the current case, this means that anomalous pairing 
gaps must necessarily be of even parity. 

We solve the system of five coupled gap equations as 
stated in Eqs. (30)-(34) numerically to obtain the fully 
self-consistent solutions in the coexistence phase. Since 
we have not included any energy dependence of the ef¬ 
fective pairing vertices, we restrict the momentum sums 
of Eqs. ([30[-(34) to states in the vicinity of the Fermi 


surface. The effective potential for pair scattering be¬ 
tween anomalous gaps, stated in Eq. (37), for¬ 

mally includes a divergent contribution for k' = k. To 
avoid such contributions we cut off all pairings r^“^’ 2 
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FIG. 8. (color online) Self-consistent solution for the normal 
and anomalous superconducting gaps as stated in Eqs. (301- 
(34 1 for the electron-doped system with /1-bands outside the 
energy cut off. The superconducting gap symmetry is d x i_ y i 
for both normal and anomalous gaps. The next-nearest neigh¬ 
bor hopping is t' = —0.35 and the bare Coulomb interaction 
is U = 3. The energy cut off set to e c = 0.15. We show the 
absolute value of the gap averaged over k-states within the 
energy range ±e c , e.g. A® a = A- |A® Q (k)|. Within this 
energy range no normal superconducting gap develops on the 


/3-band, i.e. App = 0. The dashed blue line shows A® 
case where all anomalous gaps are excluded. 


in the 


to a maximum value of y max = 200i. By this pro¬ 
cedure, we find that also the anomalous gaps acquire 
a d x 2_ y 2 structure. In Fig. [8j we present the normal 
and anomalous mean field gaps given in Eqs. (30)-(34) 
in the case of an electron-doped system in the regime 
of well-developed SDW order. We show the gap value 
averaged over k-states close to the Fermi surface, e.g. 

|A® a (k)|, where N s is the number of k- 
states with \E£\ < e c and e c = 0.15t is the cut off en¬ 
ergy. Upon approaching half filling where W = 0.85 
for this band, the self-consistent solutions gradually de¬ 
crease. This differs from the results obtained for the lin¬ 
earized gap equation in Fig. [6] where the decrease is very 
abrupt. The gradual decrease in Fig. [8] follows the de¬ 
crease in the number of states within the energy range 
around the Fermi surface towards half filling. 

From the results presented in Fig. [8j it is evident that 
the mean field A®^ is the most important contributor 
to anomalous pairing whereas A®^ and A® ^ are neg¬ 
ligible. This illustrates the point that it is the linear 
coupling to the normal gaps via the pairing interaction 
r“^ of Eq. (351 which generates the anomalous pairings. 
Without this coupling term also the subdominant mean 
fields A® A!, and A®^ would vanish. Furthermore, the 
inclusion of anomalous pairings slightly enhances the size 
of the intraband gap. This is seen from Fig. [8] by corn- 
pairing A® a as deduced from the self-consistent determi¬ 
nation of all five gaps shown by blue circles with a full 
line with the self-consistent calculation of the normal in¬ 



FIG. 9. (color online) Form factors of (a) the normal super¬ 
conducting gaps A““ and Af® and (b) the anomalous gap, 
A“A All gaps have d x 2_ y 2 structure. Whereas the normal 
gaps show the same overall phase on electron and hole pock¬ 
ets, the anomalous gaps display an internal 7r-shift between 
the electron and hole pockets. 


traband gaps alone, which is shown by the blue dashed 
line. In the latter case the gap magnitudes are slightly 
smaller. 

If hole pockets are included within the energy cut off, 
i.e. \E^\ < e c , the self-consistent calculation finds a 
7r-shift of the anomalous superconducting phase at the 
hole pockets compared to the electron pockets, see Fig. [9] 
(b). This is a robust feature arising from the momentum 
structure of the normal-anomalous pair scattering poten¬ 
tial and as a result of the internal sign between a- and 
/3-operators in Eq. (32). We have tested numerically that 
this result is insensitive to the relative ratio of normal and 
anomalous pairing potentials. 

Lastly, we mention that the self-consistent gap solu¬ 
tions are all real-valued. Therefore unlike the iron-based 
case studied in Ref. SU we do not find a TRSB supercon¬ 
ductor in the coexistence phase. In conclusion, the incor¬ 
poration of anomalous pairings in the effective interaction 
Hamiltonian gives only small quantitative changes to the 
intraband gap results presented in the previous sections. 


IV. CONCLUSIONS 

This study represents a very detailed investigation 
of spin-fluctuation-mediated superconductivity in a sys¬ 
tem with well-developed itinerant antiferromagnetic or¬ 
der and includes a calculation of the complete coexis¬ 
tence phase diagram of the Hubbard model. The spin 
fluctuations bear important fingerprints of the spin order 
with the transverse spin fluctuations corresponding to the 
Goldstone mode of the spin-symmetry-broken state. We 
find that longitudinal and transverse spin fluctuations 
are equally important for the development of an effective 
pairing glue between the fermionic quasi-particles of the 
spin-density-ordered metal. Both types of spin fluctua¬ 
tions act in concert in the even parity quasi-spin singlet 
channel and this gives rise to a robust gap solution of 
d x 2_ y 2 structure at all doping levels within the coexis- 
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tence doping region. The situation is quite different in 
the case of odd parity quasi-spin triplet superconductiv¬ 
ity. Here longitudinal spin fluctuations promote a node¬ 
less p-wave solution in the hole-doped system due to an 
effective intrapocket attraction. However, a strong in¬ 
trapocket repulsive contribution from transverse fluctu¬ 
ations destroys this effective attraction and destabilizes 
the p -wave solution. Thus, the coexistence phase treated 
within the Hubbard model does not support a nodeless 
gap on the hole doped side, in contrast to a recent study 
of the coexistence phase within the t — J model with¬ 
out double occupancy constraint.^ A modified version 
of the Hubbard model in which the relative contribution 
from the transverse fluctuations is suppressed, could also 
stabilize a nodeless p-wave state. 

Finally, we have investigated the additional interband 
pairing amplitudes that appear in the SDW phase as a 
consequence of the linear coupling to the normal intra¬ 
band gaps. Such anomalous pairing gaps are required 
to have the same parity under inversion as the normal 


intraband gaps. Specifically, the anomalous pairing am¬ 
plitudes also manifest a d x 2_ y 2 structure and the gap 
magnitude is similar to the normal pairing gaps. 
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